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Abstract 

The Roth's two-pole approximation has been used by the present authors to investigate the role of d — p 
hybridization in the superconducting properties of an extended d—p Hubbard model. Superconductivity with singlet 
d x 2_ y 2-wa,ve pairing is treated by following Beenen and Edwards formalism. In this work, the Coulomb interaction, 
the temperature and the superconductivity have been considered in the calculation of some relevant correlation 
functions present in the Roth's band shift. The behavior of the order parameter associated with temperature, 
hybridization, Coulomb interaction and the Roth's band shift effects on superconductivity are studied. 
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1 Introduction 

After almost two decade of intense research about the 
cuprates, there is still plenty of open questions in this 
problem. However, it is recognized that the electrons 
which move on the CuOi planes are the most relevant 
to describe their physical properties yQ. In the undoped 
regime, these compounds are insulators and exhibit anti- 
ferromagnetic order at sufficient low temperatures ^ |5] . 
The physical properties of the insulating phase can be well 
described by the Heisenberg model 0. Upon doping, these 
systems suppress the antiferromagnetic order and become 
superconductors. In this scenario there is no doubt that 
the d-d electron correlations play a fundamental role. 

The study of the electronic structure near the Fermi 
level Sf in such strongly correlated systems is very impor- 
tant to understand their physical properties 3 . Earlier 
angle-resolved photoemission experiments (ARPES) have 
showed the presence of flat bands close to £f in a region 
centered around the point (w, 0) in the p-type cuprates like 
Bi 2 Sr2CuOe and YBaiCu^O y 00]. Due to the presence 
of strong correlations, to study some physical properties of 
these cuprate compounds, the one-band Hubbard model 
|3] can be used. Bulut et al. E] have done Monte Carlo 
calculations in the one-band Hubbard model. Their re- 
sults show bands with a flat region near (jr, 0) point for a 
given doping which agreed with the previously mentioned 
ARPES results 0. 

Beenen and Edwards using the Roth's two-pole 
approximation 9 in the one-band Hubbard model, have 
studied the normal state of the model obtaining flat quasi- 
particle bands, which agree well with those found with 
Monte Carlo simulations [HI El- The Roth's two-pole ap- 
proximation has been proposed to improve the Hubbard- 
I approximation 5] by considering a decoupling scheme 
which produces an additional energy shift (the Roth's band 
shift) in the peaks of the spectral function. That result 
is in agreement with those obtained by Harris and Lange 
They looked at the moments of individual peaks 
in the spectral function. The presence of the exchange 
term (SiSj) in the Roth's band shift exhibits in it a spin 
dependence. As consequence, the Roth's method raises 
the possibility of magnetic solution in the Hubbard model 
while this feature is not present in the Hubbard-I approxi- 
mation. Recently, due to the good agreement between the 
Roth's and the Monte Carlo data, Beenen and Edwards 
have extended the Roth's two-pole approximation in order 
to investigate the superconducting properties of the one- 
band Hubbard model. Their main achievement has been 
to show the emergence of the pairing with d x 2_ y 2 symme- 
try in a given amount of doping. In that approach, the 
gap equation for <i-wave symmetry depends on a particu- 
lar four operator correlation function which, in principle, 
can be found extending the Roth's formalism to obtain 



two particle Green's functions. However, the authors have 
introduced two decoupling schemes to calculate the gap. 
The first one (the factorization procedure) has been for- 
mulated to treat the problem for intermediated values of 
U (the Coulomb interaction) and it provides an upper es- 
timate for the gap and T c (the critical temperature). The 
second one is adapted for very large U scenario which pre- 
serves the proper limit for U — > oo where the gap function 
vanishes. 

Nevertheless, the one-band models neglect the pres- 
ence of the oxygen sites. Due to the strong correlations at 
the Cu-sites, the oxygen sites may be occupied by holes 
when the system is doped 2 . For instance, the Hub- 
bard one-band model suffers some limitations to describe 
the low-energy physical properties of the cuprate super- 
conductors In the doped regime, the one-band Hub- 
bard model gives a wrong description of various proper- 
ties, like, for example, the asymmetric magnetic doping- 
temperature phase diagram : 12'. Therefore, a model which 
take into account also the oxygen can be more adequate 
to treat the cuprate systems in the doped regime [T3] . 
This raises the question whether it is possible to extend 
the Beenen and Edwards analysis to investigate the d- 
wave symmetry superconductivity when the hybridization 
is present. 

Recently, the present authors have used the extended 
Hubbard model ^31 with the Roth's method to study the 
role of hybridization in the superconductivity following 
closely the approach introduced by Beenen and Edwards 
|14|. As discussed in the references |3] and \H>\, the flat- 
tening of the bands is directly related to the band shift. 
The presence of flat bands at Fermi level £f in the p-type 
cuprates [3] suggests a high density of states at the Fermi 
level, which can favor pair formation. Therefore, consid- 
ering that the main responsible elements for the density of 
states at the Fermi level are the d-electrons, it has been as- 
sumed in Ref. 14 that the d—d pairs are the most relevant 
ones for superconductivity [lfij . 

The band shift plays an important role in the study 
of the superconducting properties of the model using the 
Roth's or some similar procedures. In reference |14| . the 
factorization procedure |Hj has been used to investigate 
the effects of the hybridization on the superconductivity. 
It has been shown ^1] that the hybridization has strong 
effects in the shift and, therefore, in some superconduct- 
ing physical properties such as the critical temperature 
T c . However, as a first approach, in reference [14| . the 
band shift has been evaluated taking into account the hy- 
bridization effects, but disregarding temperature effects, 
superconducting properties and, most important, it has 
been considered in the limit U — > oo. As a consequence 
of this limit, many correlation functions, which appear in 
the shift, are vanished. The important point is that these 
correlation functions are very relevant in the sense of to 
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include correctly the hybridization effects. Therefore, it 
would be necessary to calculate the shift with the U fi- 
nite in order to include the hybridization effects in a more 
complete way. 

In this work, the superconductivity problem has been 
studied using the Roth's method, following closely refer- 
ence [HI , but adapted to the d-p extended Hubbard model. 
Here, special attention is devoted to the effects of the hy- 
bridization and superconductivity in the band shift. In 
order to have the effects of the hybridization included 
properly in the superconductivity, the gap function is ob- 
tained using the factorization procedure [TI] and the shift 
is evaluated with finite U. This procedure is justified be- 
cause it preserves some correlation functions present in 
the band shift, which are non vanishing for finite U . As 
consequence, it captures the effects of the hybridization 
properly. Some preliminary results of this approach have 
been given in Ref. |17| . 

There are some shortcomings in the Beenen and Ed- 
wards approach |18l I19j . For instance, the d x 2_ y 2 pairing 
is quite dependent on the choice of the decoupling scheme 
for the correlation functions related to the gap. However, 
in the present work, the main goal is to study the effects 
of hybridization. Therefore, as discussed in the previous 
paragraph, the natural choice is the decoupling scheme for 
intermediated U, which is also the simplest one. One is 
allowed to find in that procedure, at least, a better esti- 
mate for the gap (and therefore for T c ) as a function of 
hybridization within the same decoupling procedure. 

The paper presents the following organization. In sec- 
tion it is introduced the model and given a short in- 
troduction of the Roth's method [H]. Also, some analytic 
expressions for quasi-particle bands and the Green's func- 
tions are derived. In sectional the factorization procedure 
proposed by Beenen and Edwards [S] is applied for the 
present case. In section 0] the band shift is discussed in 
detail. The numerical results are showed and discussed in 
sectional Finally, in section^ a short summary and some 
concluding remarks are given. 

2 General formulation 

The model considered here assumes overlapping bands. 
It is characterized by a narrow d-like band with a large 
density of states and a wide p-like band with low density 
of states. The extended Hubbard model is defined as: 

2,(7 hj^ i 

+ E( £ p ~ v)v\ a Vio + ]T t%v\ a Vjv 

+£#(4Pi*+pU0 (i) 



where /x is the chemical potential. The d icr> di a and pt , 
Pi a are the creation and annihilation operators of the d- 
and p-electrons, respectively, with spin a at a lattice site 
i. The £d and e p are the centers of the on site energies 
of the occupied orbitals of the copper and oxygen respec- 
tively. The second term of the Hamiltonian given in Eq. 

describes a narrow <i-band with a hopping amplitude 
t d . The Hamiltonian considers also a p-band which 
is wider than the <i-band. The following relation between 
t d and V can be established t p = at d with a > 1. Also, 
t d < to coincide the bottom of the d- and p-bands with 
the r point k x — k y — as suggested by experimental 
results |3J. The third term corresponds to the Coulomb 
interaction U that represents the repulsion between two 
holes in the same d-orbital. The last term of the Hamil- 
tonian 0) is the d — p hybridization and describes the 
nearest neighbor hopping process between the d-orbital of 
the Cu-atom and the p-orbital of the O-atom. Considering 
a rectangular two dimensional lattice, the unperturbed d- 
and p-energy bands are given by 



si = 2t d (cos(k x a) + cos(k y a)) 



and 



2t p (cos(k x a) + cos(k y a)) 



(2) 



(3) 



where a is the lattice constant. 

In this work, the Hamiltonian given in Eq. (JJJ has been 
investigated using the Roth's two-pole approximation [Hj 
to obtain the Green's function in the Zubarev's formal- 
ism. In the Roth's procedure, a set of operators {A n } is 
introduced in order to describe the relevant one particle 
excitations of the system. These operators satisfy in some 
approximation the following relation: 



(4) 



Anticommuting both sides of Eq. (@J with each oper- 
ator of the set {A n } and taking the thermal average, the 
equation (@J becomes: 



(5) 



where E nm and N nm are the energy and normalization 
matrices, given by 



[A n ,H] { _ } ,Al 



(+) 



(6) 



and 



N nm = ([A n ,Al} {+) ). 



(7) 

In matrix notation, Eq. © is written as E = K • N, 
where, if N is nonsingular, then the K matrix can be ob- 
tained. With the equation of motion (in the Zubarev's 
formalism) of the Green's function 

Gnm (w) = {{A n ;A] n )) u (8) 
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and the Eqs. J3J-0I, it is possible to obtain the following 
general Green's functions 



m ' 



(9) 



In the particular case, where B = A\ nl the elements of the 
Green's function matrix G are given by Eq. JSJ. Thus, 
using the matrices E and N, the matrix G is given by 



G (w) = G(w)N 



(10) 



where 

G(lo) = N(cjN-E)~ 1 . (11) 

Considering the fact that the operators of the set {A n } 
describe the particle excitations of the system, the choice 
of these operators is very relevant to study the physical 
properties of the system. In order to discuss superconduc- 
tivity, Beenen and Edwards, in their approach with the 
one-band Hubbard model, mixed electron and hole op- 
erators and evaluated anomalous correlation functions 8 . 
Therefore, using a set of four operators {ci a , ni_ a Ci a , cj_ CT , 
niacl_ a }, it has been obtained a four-pole approximation 
to the Green's functions. However, in order to discuss 
the role of the hybridization, it is necessary to adapt the 
formalism to include a p-operator in the original set of op- 
erators used by Beenen and Edwards. Thus, the new set 
of operators is given by 



d n a d d> n a d 



(12) 



In the present work, only the singlet pairing is consid- 
ered, and particularly the d-wave symmetry. In this par- 
ticular case, (di-a-dia-) = and '^^(di- a di <T ) = 0, where I 

i 

are the neares t ne ighbors of i. Using the set of operators 
given in Eq. Ijl2|l . ancl introducing the symmetries dis- 
cussed above, the elements of the energy matrix defined 
in Eq. © can be obtained as: 



E2 







v, 



n<L„V, 



dp 



-E 2 



Ik 



v, 



where V k p and V k a are the Fourier transform of t"J and 

t p ^ respectively. The matrix E2 present in the energy 
matrix E5 is given by: 

e d + e d k + Un d _„ (e d + e d + U)nt a ' 
(s d + e d k + U)n d _ a Un d _ a + T fe _ CT 



pd 



where Ed = £d ~~ A*- It i s assumed that the system consid- 
ered here is translationally invariant, then n_ a = n d _ a . 
The quantity T^-a is the Fourier transform of 

IV, = (s d - rintjij +4(nl CT ) 2 + n d _ a (l-n d _ a )W l3 - a . 

(15) 

In Eq. 115(1 . the band shift Wij^ a is defined as: 
_t d J {(nt a n d _ a )-{n d _ a f)+h l]ry 



WVa - " 



n d _ a {l-n d _ a ) 



(16) 



where Ay- CT can be separated into two explicit contribu- 
tions 

(17) 



A - - - \ d _|_ \P d 



The terms Af ia and A?- CT are associated with the hopping 



. V ' 

term of Ay- ff may be written as 



tfj and the hybridization t p - , respectively. The hybridized 



= E^[ 2 <P3"-X*-> - (pU*-)]^. (18) 



and the part associated to the hopping tfj is given by 
I 

+ (dl a d^ a nt a ) - {d\_J M nt a )} Sa 

~ t d 3 {{d] a d]_ a d M d^) + {d) a d\-Aj-M}. 

(19) 

The calculation of the correlation functions presented in 
Eqs. 1(18(1 and 1(19(1 will be discussed in detail in section |U 
One of the most important elements of the matrix E5 
is E24 = 7s- , where 



7fe 



mi 



ik-(Ri-Ri) 



7,7 



and 



In = (n d _ a diirdi-a + nf a d la di^ a ). 



(20) 



(21) 



The correlation function j k gives the gap of the supercon- 
ductor state in the d— wave case. 

The elements of the normalization matrix N5 are given 
from Eq. J7J as: 



(13) 



N n = N 33 = N 55 = 1 



.dp 



and 



N12 = N 21 = N 22 = N u = N 43 = N 4i = n°_ 



(22) 



(23) 



E/9 



(14) 



The remaining elements of the normalization matrix N5, 
due to the rf-wave symmetry and the anticommution rules, 
have been found to be zero. 



4 



with 



Using the energy and the normalization matrices E5 where 
and N5, respectively, the matrix Green's function G5 de- 
fined in Eq. I|10|) can be obtained. For simplicity, only the 
most relevant elements (for the purposes of this work) of 
this (5 x 5) G5 matrix are shown. Following the Roth's 
notation [S], the correlation function (BA) is related to 
the Green's function ((A;B)) U as: 



£H = 2?H-7 fe V-£i 2 i) 



(34) 



(BA) = T U ((A;B)) U = ±-l duf(u){(A; B)) u> (24) 

where /(w) is the Fermi function. The chemical potential 
/i is obtained in the standard way, using the clement Gi^ 
of the matrix G5 and the relation given in Eq. Ij24|l . The written as: 
matrix element G- 1 is given by 

ka 



V(lo) = [{u - E n ){un d _ a - E 22 ) - {umtr - E 12 ) 2 } 

x [(lu + Eu)(u}ni a + E22) - {un d _ a + E 12 ) 2 ]. 

(35) 



In Eq. (|35(l . E\ 2 and £"22 are elements of the energy matrix 
E2 given in Eq. (|14|) . The use of a set of five operators A n 
results in a five-pole approximation to the Green's func- 
tions. Then, the D(uS) defined in Eq. (|33|l may be also 



[u~E bb )[A{uj)-{uj + E ll )^ k 2 



(25) 



D(iu) = (n d _ a ) 2 (l-n d _ a ) 2 {LJ-E lk )(u-E 2k )(u;-E 3k ) 
x(uj- E 4k )(tu- E 5k ) (36) 



where the quasi-particle bands E pk (with p = 1, .., 5) sat- 
where En and E 55 are elements of the energy matrix E 5 , isfy D = det(ujN 5 - E 5 ) = 0. Therefore, the resulting 
defined in Eq. I^J- In Eq. it is also necessary to Green's function can be written as a sum of five terms: 

introduce the following definitions: 



A(w) = (nl a ) 2 (l-ni CT ) s 



G£M=£-^ 



xH+4V+«r>+<o (26) 



(2) 



S 3 )\ 



E ,- 

pk<r 



(37) 



where Z ^ a express the spectral weights which satisfy 



with 



„( 2 ) 



and 



7 (2) 
^ka ^ka 



— En, 



(2) 



En 



a 



(3) 
ka 



-z^zf\z^ 

ka ka \ ka 



z(2) 
J ka 



En). 



The quantities Z^J and Z y k J are defined as 



(2) 



(i) = U + 2(e d -fi) + e d k + W k - a 

ka o 



Ha 



and, 



7 (2) _ o-(l) 

ka — ^ka 



Ha- 



(27) 
(28) 

(29) 



(30) 



(31) 



lka 



Z 2ka 



Z 3ka 



Z 4ka 



Z K t = 1. 

oka 



(38) 



3 Calculation of the gap function 
using the factorization procedure 

In the case of d-wave symmetry, the traditional correlation 
function (di^^di^) is always zero. Therefore, this correla- 
tion function can not be used to determinate the pairing 
gap in the d-wave channel jHl^j- In the factorization pro- 
cedure proposed by Beenen and Edwards in Ref. [B], the 
correlation function given by Eq. (|21J) is rewritten as 



In the particular case, when Ed, j k and t pd are zero, Z^J 
(2) 

and Z ka represent the quasi-particle bands in the para- 
magnetic normal state of the one-band Hubbard model. 
The term is given by: 

&ka = \J(U + W k - a - e d f + 4n d _ a U(e d - W k ^ a ) (32) 

where W k - a is the Fourier transform of Wij- a given in 
Eq. (|16() . The denominator of the Green's function Gi 1 



Iti = [(4-a d i-<?) + (4 l7 d i<J )}(d l - lT di a ) 



(39) 



where the symmetry 7 i; = 7 ;i is conserved, and the prod- 
dirr are split up. It is also introduced 



ucts d\ a di- a and d, 



n. 



i — a^ia 

oi. = (4- 



which allows to rewrite Eq. I|39|l as 

2n d i a (d i ^ rT di, 



given in Eq. I|25|) is defined as: 

D{u>) = (u;-E 55 )D(u,)-V k dp V k pd [A(u;)-(u, + Eii)j k 2 l 

(33) 



la 



where rim„ can be calculated from G\\ 



(40) 



(41) 



01o . 0011 ^a^uiarcu iium v_r fe(J . 

Considering the c?-wave symmetry, the Fourier trans- 
brm of 7 iZ given by Eq. (121)1) becomes 



Ik = 9 [ cos (kx®) — cos (k y a)] 



(42) 
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where 

9 = 2^7 (43) 

is the gap-function amplitude. Due to the <i-wave symme- 
try, 7 i; = +7 for Ri — Ri in the x direction and 7^ = —7 
when Ri — Ri is in the y direction. The Fourier transform 
of the correlation function (di- a di a ) is given by 

(di-^dia) =^J2 <? HRl - Ri) (dk-M (44) 
k 

where L is the number of sites in the system. The cor- 
relation function (dk-a-dka) can be evaluated using the 
Green's function and the relation given by Eq. J21}. 
The Green's function GH can be rewritten as: 



where 



G%(w) = -%U*Fg (45) 
(n d _ a ni-n d _ a ) 2 (u;-E 55 ) 



D(u) 



(46) 



and D(lo) is defined in Eq. i|33[l . 

Combining the equation (|2()|l with the Eqs. (|41|l to 
l(45[l . the gap equation can be written as: 



i k = -iM la t d u 2 i a 



where 



(47) 
(48) 



with 



h = -^-<j>f(u)F a (u})duj 



4 Definition and calculation of the 
band shifts 

Using the definition (|17|) in Eq. JTBJ, the band shift Wij a 
can be written as: 



Wu n = W d + W pd 



(50) 



where 



tfMnf „n d „) - (n d <T ) 2 \ + Af,„ 



l]-<T 



n d _ a (l-n d _ a ) 



and 



W, 



pd 



A pd 



n d _M - n d _ a ) 



(52) 



The quantity A pd a is given by Eq. fH^I . The correlation 
function (pj_ a di- a ) present in A?. ff can be obtained from 
de Green's function 



[A(co)-(u J + E 11 )j k 2 ]V pd 
D(cj) 



(53) 



The remaining correlation function (pj_ nf^di—o) present 



in A P f„ is calculated from the Green's function 



where 



nU [B(w) ->+£ ; u )7 fc 2 ] V dp 
D(u) 



(54) 



B (lo) — A (lo) + n_ CT (l — n_ a ) UT>i{lo) (55) 

with A(u>) defined in Eq. l|26l) . The quantity X>i(o>), in 
terms of the elements of the energy matrix is given 
by: 

2?i(w) = (w - £u)(aml ff - £ 22 ) - (wnl, - £ 12 ) 2 . (56) 

The Green's function G 2 ,^. tends to zero as U — > oo, 
consequently, the correlation function (p\_ a n d a di- a ) also 
vanishes recovering the result of Ref . for A^ rf CT . 

The quantity Af J(T present in Eq. lj5T|) is given by Eq. 
H19|) . The Fourier transform of Wf^ a is given by: 



(57) 



Substituting Eq. 1)19)1 into Eq. 1)51)1 and then putting the 
result into Eq. 1)57(1 ■ the Fourier transform of Wh_ a can 
be written as: 



X<e ifc '^{{n^nt)-(nt) 2 



+ 



(4 <r d J -_ <r 4_ (T do ff ) - (4a4-^*-" d <"'>}(' 58 ) 



The correlation functions present in W d a are evaluated 
following the original Roth's procedure jH). Introducing 
extra operators Bi a , the correlation functions of the form 
(A n B i(7 ) can be calculated by using Eqs. © and In 
Refs. jHlEl) the sum present in Eq. © has been con- 
sidered only over the operators which describe the normal 
state of the system. In the present work, the sum includes 
also the hole operators which describe the superconduct- 
ing properties of the system. Thus, W d a is given by: 

nj(l nj)Wl = h la + Y. *o/^ (H* + h V°) ( 59 ) 

where the term h^ju is directly related to the gap function 
7r through the Green's functions Gi 3 and G- 4 (see Ap- 
pendix A). The quantities ft-io-, /i 2 jo- and ft^cr are given in 
Appendix A. 
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Table 1: The different approaches considered to evaluate 
the band shift WV . 

Acer 



Beenen and 
Edwards [E] 

Ref.[Hj 

Present work 



5 Results 

In this section, the numerical results obtained in this work 
are presented. One of the most important parameters of 
the model given in Eq. (JTJ is the d — p hybridization 
|14l |23] , which is defined as 

V~ p = -iV dp [sm(k x a) ~ sm(k y a)}. (60) 

In this work, as in Ref. |20|. the hybridization has been 
assumed fc-independent (V dp ) 2 = (V~ p VS d ), where (...) is 
the average over the Brillouin zone. In Ref. [HI Sengupta 
and Ghatak have also used a fc-independent hybridization 
due to the fact that the pairs occur within a small energy 
interval around the Fermi level, therefore the dispersion of 
the hybridization can be neglected. 

The total occupation number is given by nx — n d + 
n d _ a , where n d is obtained combining G^. (Eq. lf^5|l ) and 
the relation given in Eq. Ij24(l . The charge transfer energy 
A = e p — Ed is positive. This means that the first hole 
added to the system will energetically prefer to occupy 
the (i-orbital of the copper ions p^. All results presented 
in this section are obtained with Ed = and e p = 3.6eV. 
Consequently, A = 3.6eV, as estimated in Ref. |22| . 

As discussed in Refs. 010], the band shift (see 
Eq. (JSUJl) can be evaluated considering different approx- 
imations. In the limit U — * oo, some terms of the band 
shift vanish (see Ref. (Hj). In Ref. the present au- 

thors estimate Wr in the limit of U — > oo, but with finite 

KG > 

U in other parts of the problem. In Ref. [Hj, Beenen 
and Edwards evaluated W% in the normal state (where 
7^: = 0) and considering T equal to zero and finite U us- 
ing the one-band Hubbard model (hybridization null). In 
the present work, the correlation functions present in Wt 
given in Eq. (|58|) are evaluated following closely the pro- 
cedure used by Roth in Ref. Nevertheless, here, also 
the hole operators given in the set of Eq. l(T2"|> are used 
to evaluate the correlation functions. As consequence, a 
new term (h S j a ) appears in (see Eq. (JSSJ). The 

approximations used to evaluate are shown in Ta- 

ble H In figure da), the quasi-particle bands E p k, with 
p = 1..5 (see Eq. (jbSIl '). are plotted along the symmetry 




k k 



Figure 1: (a) The quasi-particle bands for U — I2\t d \, 
Vg d = 0.2\t d \ and n T = 0.76. (b) The electron and hole 
bands in the neighborhood of the (tt, 0) point where the 
gap structure is relevant. The dashed lines show the result 
for V pd = 0.0 while the solid bf lines correspond to Vq = 
0.2\t d \. (c) shows the hybridization (V£ d ) gap near the 
point (7T, 7r). 

lines (0, 0) — (tt, it) — (tt, 0) — (0, 0), in the two-dimensional 
Brillouin zone. The quasi-particle energies E p k, in the su- 
perconducting state, are relative to the chemical potential 
[i. The circles bf show the e\ band, where the center of 
is shifted by e v = 3.6eV relative to the zero of energy. 
All results shown in this paper are obtained with t p = 2t d . 
The dashed line corresponds to the noninteracting (U = 0) 
band si relative to the noninteracting chemical potential. 
The figure Gib) shows the superconducting gap between 
the electron and hole bands in the neighborhood of the 
(tt, 0) point, while on the k x = k y diagonal (Fig. Wi a )) the 
gap is zero. This fact reflects the d-wave symmetry pro- 
posed in this work. The dashed lines show the absence of 
the gap in the normal state. In figuredc), the region near 
to the (tt, 7r) point shows the gaps produced by the d — p 
hybridization Vjj pd . The dashed lines show the result for 
Vq 4 = 0. In figure |2 the electron and hole quasi-particle 
bands are shown for two different values of hybridization. 
As can be observed, the hybridization shifts the quasi- 
particle bands to lower energy by breaking the symmetry 
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Figure 2: The electron and hole bands in the region close 
to the (tt,0) point for U = 12\t d \ and n T = 0.76. The 



solid lines correspond to Vq 4 = 0.0, while the dashed lines 
show the result for K pd = 0.3\t d \. 



in relation to the k axis. 

The figure 01 shows the spectral weights Z % a for two 
different hybridizations. The dashed line corresponds to 
the sum of the five spectral weights which is equal to one 
(see Eq. In figure|2Ib), the effects of the hybridiza- 

tion on the spectral weights are shown. Such effects cause 
a small change in the chemical potential and consequently 
in the superconductivity. 

The figure U shows the behavior of gap function ampli- 
tude j as a function of the hybridization VQ d . It is clear 
that there is a decreasing of g with increasing V^ pd . 

The analysis of the function F a {ui) introduced in Eq. 
(|48|l and defined in Eq. (|49|l is important to understand 
the behavior of the gap function amplitude showed in fig- 
ure 01 The figure shows the function F a {ui) for T = 
and two different values of hybridization. As can be seen 
in the dashed line, the magnitude of the function F a (u>) 
decreases when the hybridization is enhanced. Moreover, 
the function is shifted to lower energy, breaking the sym- 
metry respect to lo = 0. The symmetry break has been 
also observed in figure El for the electron and hole quasi- 
particle bands. For T = 0, the product f{w)F„(tJ) given in 
Eq. <|48(l vanishes when to > 0. That is because the Fermi 
function f(uj) is zero for that range of u>. As consequence 
of the shift and the suppression of F a (uj), the value of I a , 
which is given by the integral in Eq. H48|l , decreases when 
the hybridization increase. However, from Eq. (|47|) . it is 
necessary a minimum value for I a to obtain a nonzero so- 
lution for 7. But, for very strong values of hybridization, 
the minimum value for I a is not reached and only the zero 
solution exists. 

According to this analysis, there is a critical value of 




(0.0) 



(0,0) 



Figure 3: The spectral weights Z^ a for U = 12|i d |, tit = 
0.76, T = and two different values of hybridization. 

hybridization (Vq c ), above which, the superconductivity 
is suppressed. Similar results, which show a critical value 
for the hybridization, were also obtained in Ref. [221, for a 
fc-dependent hybridization and using the Hartree-Fock ap- 
proximation for the electron-electron interaction. In Refs. 
|16l I2.'ij . although the high T c was not considered, the 
hybridization effects play an important role for resonant 
states. The discussion above is also valid if the values of 
the temperature T are raised with V pd constant. The only 
difference is that in this case the Fermi function becomes 
sloping smoothly, changing the product f(uj)F a {uj). The 
effect of the temperature in the Fermi function causes a 
decreasing of l a and consequently of T c . 

Since the hybridization is directly related to the ap- 
plied pressure , the transition temperature T c may have 
a dependence on pressure through the hybridization. How- 
ever, the pressure dependence of T c is very complicated in 
high temperature superconductors. As discussed in Ref. 
|16| . at least, in conventional superconductivity where the 
electron pairing is mediated by phonons, two effects are re- 
sponsible for the pressure dependence of T c . The first one 
is related to the lattice vibrations, while the second one 
comes from the electronic contribution. As long the pres- 
sure is increased, the lattice vibrations tend to increase 
T c , whereas the effects of the electronic contribution as- 
sociated with the hybridization cause a decreases of T c . 
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Figure 4: Behavior of the gap function amplitude g as a 
function of the hybridization V^f d for T=0.004eV. 
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Figure 5: Function F a (u) (defined in Eq. I|49|)'l for U = 
12|i d |, tit = 0.76, T — and two different hybridizations. 

The figure shows the function amplitude g, as a func- 
tion of temperature T for several values of hybridization. 
If the hybridization is enhanced, the critical temperature 
T c decreases. Therefore, the present results agree with 
the discussion above in the scenario where the electronic 
effects dominate. This behavior for T c is also shown in 
figured in phase diagrams displaying T c versus the total 
occupation number nr. The numerical results obtained 
show that there is a critical value of hybridization where 
the superconducting phase vanishes. These results agree 
with the ones obtained by the authors of Ref. 24 for 
heavy-fermion superconductivity with an X-boson treat- 
ment. 

The solid lines in figures Ufa)- (b) show the present re- 
sult for T c , where the effects of the temperature, the su- 
perconductivity and the Coulomb interaction have been 
included in the calculation of the band shift Wt ■ The 
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Figure 6: The gap function amplitude g, as a function 
of temperature T, for U — 12|< rf |, jit = 0.76 and several 
values of hybridization. (t d = — 0.5eV.) 



dashed lines correspond to the results obtained in Ref. 
|14| . where the band shift has been evaluated consider- 
ing T = 0, 7£ = and U -> oo (see table |TJ). The 
difference between the results can be explained by the 
analysis of the Eq. I|18|l where some of the correlation 
functions present in the band shift vanish in the U — > oo 
limit. In Eq. . the correlation function (pj_ cr nf a .di- a ), 
which is directly related to the hybridization effects in the 
band shift, vanishes for U — > oo. It is important to high- 
light that the correlation functions in Eq. 1)18(1 are both 
negative. Therefore, for large U, the correlation function 



(Pi 



decreases and the hybridized shift Wj? is 



enhanced. However, for intermediate values of U, both 
correlation functions remain finite. As consequence, the 
hybridization effects in the band shift and, therefore, in 
the superconductivity, are weakened. The figures 0c)-(d) 
show the present results when the value of U is increased. 
The main consequence is, within the factorization proce- 
dure, to shift the window of doping where superconduc- 
tivity is found, as in Ref [Sj. 

In figurelHIa), the chemical potential is show as a func- 
tion of the total occupation number tit for U = 12\t d \ and 
two different hybridizations . In Ref. [21], the authors 
criticized the Roth's method because the compressibility 
k = ^jjf- is negative in the vicinity of half-filling in the 
Beenen and Edwards result. In Ref. ^Sli by using a com- 
posite operator approach and imposing the Pauli principle, 
the authors have showed that the compressibility remains 
negative. However, they also showed that the pairing de- 
creases the strength of the negative compressibility. 

In the present work, a careful study about the nature 
of the negative compressibility and the effect of the hy- 
bridization near half-filling in Roth's approximation has 
been carried out. It has been verified that the most im- 
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Figure 7: T c as a function of the total occupation number 
jit- In (a) and (b), the dotted lines show the previous 
results from Ref. [H] for U = 12\t d \. The solid lines 
show the behavior of T c in the present approach. The 
figures (c) and (d) show the present results for U = 8\t d \. 
{t d = -Q.heV.) 



portant contribution to provide negative compressibility 
comes from the spin-term (SjSi) present in the c?-part 
Wi of the Roth's band shift W £<7 (see Eqs. El and 
IA.34(I . In reference [H], it has been showed that the cor- 
relation function (SjSi) plays an important role on the 
flattening of the quasi-particle bands. The correlation 
function (SjSi) increases with occupation and its effect is 
pronounced near half- filling. Nevertheless, when the hy- 
bridization is present, the numerical results show that it 
acts in the sense of suppressing the negative compressibil- 
ity near half-filling. Because the hybridization considered 
here is fc-independent [20], the hybridization term W£ d of 
the band shift is constant within the Brillouin zone. Its 
main effect is to shift the poles of the Green's functions 
and consequently to change the value of the chemical po- 
tential suppressing the negative compressibility. In figure 
IBIa), it is clear that the effect of the hybridization in the 
chemical potential decreases the negative compressibility. 

The figure |5fb) shows the gap function amplitude ~g 
as a function of the total occupation number. This result 
agrees with those obtained in figure ^ where g decreases 
with increasing of Vq . 

6 Conclusions 

In this work, the Roth's two-pole approximation is ex- 
tended to study the superconducting properties of the ex- 
tended Hubbard model given in Eq. Q). The quality of 
the Roth's two-pole approximation had been investigated 
in a previous work by Beenen and Edwards 8 . In their 
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Figure 8: (a) The chemical potential as a function of the 
total occupation number for U = I2\t d \, kT = 0.004eV 
and different values of hybridization, (b) The gap function 
amplitude, as a function of the total occupation number, 
with U and kT identical to (a). 



work, they showed the remarkable agreement between the 
Roth's and the Monte Carlo results for the one- 
band Hubbard model in the paramagnetic normal state. 
Moreover, the flat bands obtained with Roth's procedure 
show a qualitative agreement with the ARPES experiment 
data P] in cuprates. It is important to point out that the 
flattening observed in the quasi-particle bands which pro- 
duces a peak in the density of states, can be connected 
with the Van Hove scenario. In cuprate systems the Van 
Hove singularity is present in the vicinity of the Fermi en- 
ergy. Therefore, it is believed that the Van Hove scenario 
play a fundamental role in order to clarify the mechanism 
which drives the transition to superconductivity in these 
interesting materials |25| . 

The accuracy of the Roth's results is very related to the 
adequate evaluation of the band shift. Therefore, the focus 
of the present work has been to evaluate the Roth's band 
shift taking into account relevant effects as Coulomb inter- 
action, temperature, superconductivity and hybridization. 
Also, the effect of the hybridization in the superconducting 
of the model has been studied. This work has been carried 
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out following the factorization procedure proposed by Bee- 
nen and Edwards [S] . In order to study superconductivity, 
Beenen and Edwards proposed to include hole operators 
in the original set of operators that describes the normal 
state of the system. These operators can introduce the 
pairing formation in the d-band. The factorization proce- 
dure proposed by Beenen and Edwards [S] and the d-wave 
symmetry are considered to obtain the gap function am- 
plitude. The hybridization effects are considered by also 
including a p-operator. Thus, the set of operators is en- 
larged to five, which results in a five-pole approximation 
to the Green's functions. 

The hybridization effects present in the band shift come 
from some correlation functions. The important point is 
that part of them vanish when U — > oo, as it have been 
done in Ref. 11. In order to consider properly the hy- 
bridization effects, the band shift should be obtained for 
finite U. In fact, the obtained phase diagrams show that 
the presence of superconducting order exists in a larger 
range of doping when compared with the U — > oo limit 
|14| . for the same hybridization. Therefore, this result 
suggests that, in the U — > oo limit, the hybridization ef- 
fects are overestimated. That is the ultimate justification 
for the use of the factorization procedure :8. , which is valid 
for intermediated values of U for the gap function. 

The Beenen and Edwards's [H] results are recovered 
taking V$ d — in the present work. The hybridization 
VQ d breaks the symmetry between the electron and hole 
quasi-particle bands, respect to k axis. Also, the gap am- 
plitude function g and the critical temperature T c are sup- 
pressed with increasing the hybridization V^ pd . The results 
show that the chemical potential does not change signifi- 
cantly away the half-filling. However, near half-filling, it 
is showed that the negative compressibility decreases with 
increasing Vq . The correlation functions present in the d- 
part of the band shift Wt were discussed in detail. When 
the hole operators are also considered to obtain this cor- 
relation functions, a new term appears in the d-part of 
the band shift Wt ■ The new term is directly associated 
with the superconducting properties of the system. Nev- 
ertheless, this term is quite small and therefore may be 
disregarding in the calculation of the band shift. 
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A Appendix 

The correlation functions present in the band shift Wi 

kcr 

can be evaluated by introducing extra B operators, as in 
the original Roth's procedure. Combining the Eq. JjJJ and 
the relation given in Eq. I)24[l. it is possible to write 

(BA n ) = Tuj ^2 Gtim(w)([A m ,B]( + )), (A.l) 

m 

where A n and A m are members of the set of operators 
given in Eq. |T^|. For evaluate {nf^n^) — (n^) 2 , it has 
been necessary to introduce the following B operators: 

*£ = -^E«-*MiX (A.2) 

and 

*gl = 7^£^*<i-,4,. (A.3) 

i 

By considering the operator given in the Eq. 1A.2|) . 
the correlation function (n^n^) can be written as: 

k 

where the right side of Eq. (|A.4Jl may be obtained using 
the relation given by Eq. I)A.1|) . Therefore, it is necessary 
to evaluate the anticommutators [A m , ]( + ) for the set 
of operators A m given in Eq. (|12|l . For m = 1..5, the A 
operators are given by: 

A^ = ^£e*V-,An (A.6) 

and 

Thus, the following results have been obtained 

([^ife. = n °° ~ ^ < rf o<^v>, (A.10) 

(ATI) 

<[^4i](+)>= ' (A. 12) 
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a^> *£!,](+)> = -K4A-*), (A.i3) 

where, it has been assumed that the brackets are real and 
unchanged when the indices and j are interchanged. 
Also, due to translational invariance of the system, n,Q a = 
n d a . Considering the relations given by Eqs. (|A.1|> and 
1A.4I) with the results from Eq. (|A.10I) to Eq. (|A.14|) . the 
correlation function (n d a nQ a ) can be written as: 



ja u 0-a a 0a/ 



(A.15) 



where n?. 



the following definitions: 



i d . In Eq. 1A.15II . it has been introduced 



l 0ja — {docrdja) — — ^^J 7 , M\ : 



11 e ifc--Rj 
ka 



12 ifc-it,- 



- 7 51 T " G \ 



11 gife-Rj 



12 gifc-fij 



and 



14 gifc-itj 



(A.16) 

(A.17) 
(A.18) 
(A.19) 

(A.20) 



where Gi 1 is given in Eq. (|25(l . The remaining Green's 
functions Gi 2 and Gi 4 are given respectively by 



and 



n d _ a (uj - E 55 ) [Bjcj) - (lu + ff n )7 fc 2 
D{u) 



(A.21) 



G^H = {n d _ a f{l-n d _ a fU lk 

{u - E 55 )(u + En - Urf 



• (A.22) 



D(lo) 

where B (uj) is defined in Eq. l|5f)|) and D (ui) in Eq. (jSEJ. 

It is also necessary to define 



Gi 1 M - G 1 - 2 M 

P^ll / x _ ka K > ka K > 

GzAu) - 



l-n d 



G l 2 M ~ n -<r G V H 

Gi 2 (w)= fcCT 



and 



G tM 



n d _ a {l-n d _ a ) 
Gl^)-n d _ a GfJco) 



n d _ a {l-n d _ CT ) 



(A.23) 
(A.24) 

(A.25) 



where Gi 3 is given in Eq. 145|) . 



The correlation function 



L t)-a IL ja 



d ja ) present in Eq. (|A.15j) . 



can be obtained by repeating the procedure above using 
the operator i?l 2 . (given by Eq. Thus, 

= ot a n d _ a + p a (n d _ a n d _ a ) + af}n$ a 



where 



> {1) 



13 e ik-Rj 
ka 



14 Ak-R-j 



and 



with 



13 - - \^ t a 1 '- 



13 e ik-Rj 
ka 



Gi 3 H-Gi 4 (o;) 



(A.27) 
(A.28) 

(A.29) 
(A.30) 



Reversing the spin labels i.e., a — > — a in Eq. (|A.26|I 
and substituting the result into Eq. I|A.15(I . then 

d d OL a n a — a.j a riQj a + P a OL^ a n a — f3j a rrij a 

\ n ja n 0a) — 



1 



1 - PaP 

(!) J 1 ) 



(A.31) 



1 - PaP-a 
HPaP^l + 0^){nU\-Aa) 



For evaluate the two last correlation functions present 
in Eq. I|58|) , the following operators have been introduced 



B^ 



kja 



-ik-R, A 



i+ja u i+3-v u i-a 



and 



i-\-jcr i-\-j — tr * f" ' 



(A.32) 



(A.33) 



Using i?3, and following the procedure outlined above, the 
correlation function (dj t7 dj- <7 d _ cr dQ tT } is given by 



(SjS ) = (d\ a dj- a dl_ a d 0a ) 



l + Pa 



[a j(7 7 



Oj-a 



(A.34) 



Similarly, using B4 



(d\ <J d]_ a d Q - IJ doa} 



+ 



I -Pa 
fc ] «ad]gdU) 
I -Pa 



(A.35) 
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where the d-wave symmetry has been considered, there- 

fore, (d]J 3 _ a ) = 0. h ^ 

The four B<» operators introduced up to now are ex- 
actly the same operators used by Roth in Ref. [5] to ob- 
tain the band shift Wk* in the normal state and without 
hybridization. However, in the present work, due to the 
presence of the hole operators (see Eq. (fT2*|L a new B 
operator, which is given by 



l 3j* 



B (6) = 1 y 



ik-R z J 



diudi+jadi+j-,! , 



(A.36) 



has been introduced. With this operator, the correlation 
function (ng .d] .c?]_ .) present in Eq. (|A.35|) may be eval- 
uated. Thus, 



with 



(doo-dj-adjcrdocr) — 



J* 



\ 1 - 0*0-* 1 + 0* 

a J* n Qj-* + 0J*(j4j-* ~ m J-<y) \ 

1-0* J ' 

1 - /?«r/?_«r 
(1) (1) , (1) 

1 + A, 



(A.42) 



"jo- 



l-0 a 1 - 



(A.43) 



(A.44) 



(A.37) 



Substituting the result (|A.37(1 into Eq. I|A.35(1 . the corre- 
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{d] a d]_ a do- a d 0a ) 
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1-0* 

{ * ] \4* ] <* + & «* - [3: 

[4 



(A.38 



The result given in Eq. (|A.37(1 can be used in Eq. 
ljA.31|) to obtain 

/n d n d \ - (n d \ 2 a i an °3° + $0°™]° 
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